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TRIANGULAR BASES OF INTEGRAL CLOSURES 


HAYDEN D. STAINSBY 


Abstract. In this work, we consider the problem of computing triangular 
bases of integral closures of one-dimensional local rings. 

Let ( K , v) be a discrete valued field with valuation ring O and let m be the 
maximal ideal. We take f e 0[x], a monic irreducible polynomial of degree n 
and consider the extension L = K\x\/(f(x)) as well as Ol the integral closure 
of O in L, which we suppose to be finitely generated as an O-module. 

The algorithm MaxMin, presented in this paper, computes triangular bases 
of fractional ideals of Ol- The theoretical complexity is equivalent to current 
state of the art methods and in practice is almost always faster. It is also con¬ 
siderably faster than the routines found in standard computer algebra systems, 
excepting some cases involving very small field extensions. 


1. Introduction 

Let ( K , v) be a discrete valued field with valuation ring 0. Let m be the maximal 
ideal, irema generator of m and F = 0/m the residue class field. 

Let K v be the completion of I \, and retain v : K v —> Q the canonical extension 
of v to a fixed algebraic closure of K v . Let 0„ be the valuation ring of K v . 

Let / 6 0[x ] be a monic, irreducible polynomial of degree n and fix a root 9 e K 
in an algebraic closure of K. Let L = K{6) be the corresponding finite extension 
of K and let Ol be the integral closure of O in L, which is a Dedekind domain. 
We denote the set of non-zero prime ideals of Ol by V. 

We suppose that Ol is finitely generated as an 0-module. This condition holds 
under very natural assumptions; for instance, if L/K is separable, or ( K , v) is 
complete, or O is a finitely generated algebra over a field ■T21 Ch.I, §4]. Under this 
hypothesis, Ol is a free 0-module of rank n. An 0-basis of Ol is called a v-integral 
basis of Ol- 

The computation of v-integral bases is an essential task in computational arith¬ 
metic geometry. We are interested in constructing reduced triangular v-integral 
bases. Let us clarify these concepts. 

Definition 1.1. For each prime ideal p e V, we consider the normalised valuation: 

w p := e(p/m) -1 Vp : L —> e(p/m) -1 Z u {oo} , 

where v p is the canonical discrete valuation of L attached to p and e(p/m) is the 
ramification index. Also, we define 

w: L — > Q u {oo} , a Min {w p (a): p e V} . 
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Clearly, an element a £ L belongs to Ol if and only if w(a) ^ 0. 
Definition 1.2. A triangular family of elements in Ol, are elements 


9o{0) 9i{0) 9n-1(0) 



such that for all 0 < i < n the polynomial gi(x) £ 0[x ] is monic of degree i, and ki 
is a non-negative integer such that ki ^ w ( gi(0 )). 

A triangular basis of Ol is a triangular family which is a v-integral basis. 

Definition 1.3. A family ai, ..., a n £ Ol is called reduced if for any family 
Q-i ; ... ; a n £ O. 



There are well-known conditions on a triangular family characterising when it is 
a reduced basis. 

Theorem 1.4. Let B = {go{0)/n^ Vo f ..., g n -i{9)/'K^ Vn - 1 ^ be a triangular family 
in Ol with z/j = w{gi{9)) for all 0 < i < n. Then, 

(1) B is a triangular basis if and only if [zyj > [re (h (0))J for all h £ 0[x ] monic 
of degree i and all 0 < i < n. 

(2) B is a reduced triangular basis if and only if Vi > w ( h(9 )) for all h £ 0\x\ 
monic of degree i and all 0 ^ i < n. 

Therefore, to construct a reduced triangular 0-basis of Ol it suffices to find for 
each O^Kna polynomial gi £ 0\x\ monic of degree i such that w (gi(6)) is 
maximal amongst all monic polynomials in 0\x\ of the same degree. 

The aim of the paper is to present the MaxMin algorithm, a very simple procedure 
to construct these optimal polynomials gi. 

The desire to compute triangular local bases comes from their utility in con¬ 
structing global bases. Let A be a PID and let B be the integral closure of A in a 
finite extension of its field of fractions. If B is free as an T-module, then an T-basis 
of B can be computed by patching local integral bases of B m as an T m -module for 
an adequate finite set of non-zero prime ideals m £ Spec(T). 

This patching process, usually based on the Chinese remainder theorem, is effi¬ 
cient only if the local bases are triangular. In this case, the global T-basis of B one 
obtains is triangular too. 

The property of reducedness is important for certain applications in function 
fields .1]. The MaxMin algorithm has the advantage of producing local bases which 
are also reduced. 

MaxMin is an OM algorithm, it works with data derived from an OM factorisa¬ 
tion of the polynomial / in O v [x\. In 1999, J. Montes extended some ideas of Ore 
and MacLane and implemented an algorithm to compute a representation of the 
prime ideals of Ol by way of factoring the defining polynomial / over O v \x\. This 
“Montes algorithm” coupled with work by K. Okutsu on constructing explicit in¬ 
tegral bases of local fields, gave rise to several theoretical developments concerning 
OM representations of prime ideals [ 3 ] m m m m m mj. 

There are other methods for the computation of integral bases based on a pre¬ 
vious computation of an OM factorisation of the defining polynomial /. In [6], Sec. 
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6 ] a method was presented based in the computation of certain multipliers , follow¬ 
ing an old idea of Ore. In [7] a more direct method of the quotients was presented, 
which obtains the numerators of a reduced basis by multiplying certain polynomials 
obtained as a by-product of the OM factorisation algorithm. 

These OM methods are extremely fast in practice and their theoretical com¬ 
plexity is lower than that of the traditional methods based mainly on the Round-2 
and Round-4 routines by Zassenlraus and Ford. However, both OM methods yield 
non-triangular bases, and so a triangularisation routine must be applied to the local 
bases before they can be used to construct a global basis. These linear procedures 
are slow in practice and constitute a bottleneck for the whole process. 

The MaxMin algorithm yields reduced triangular local bases by a direct method, 
which avoids the use of linear techniques. This makes MaxMin much more efficient 
in practice. Another advantage of MaxMin is that the method is equally valid 
for the computation of bases of fractional ideals. This is particularly useful for 
the computation of bases of Riemann-Roch spaces of divisors of algebraic curves, 
which requires the computation of bases of certain fractional ideals attached to the 
divisor. 

It is common in many computational algebra systems, to provide bases in Her- 
mite Normal Form (HNF). This serves two purposes, the first is that canonical 
bases simplify the comparison of the rings that they generate. The second is that 
bases in HNF are triangular, and so patching of global bases from local bases is 
more efficient. However, the routines used to compute the HNF of a given basis 
require considerably more time than the computation of a basis. 

By using the MaxMin algorithm we can offer two distinct improvements. By 
computing a triangular basis directly, in many circumstances, we do not need HNF 
at all, resulting in a significant improvement in execution time. Secondly, if HNF 
is indeed required, it is faster to compute the HNF of a basis which is already 
triangular, compared to a random basis. 

We review OM representations of prime ideals in Section [2] In Section [3] the 
algorithm MaxMin will be introduced. In Section [|] we will discuss some computa¬ 
tional examples. Finally, the proofs of the two main theorems from Sections [2] and 
[3] will be deferred until Sections [5] and [6] 

2. OM REPRESENTATIONS OF PRIME IDEALS 

The prime ideals of Ol are in 1-to-l correspondence with the prime factors of 
/ in O v [x\. Let / = n pe -p Fp be the factorisation of / into a product of monic 
irreducible polynomials F p e O v \x\. Let n p = degF p = e(p/m)/(p/m). 

Inspired by ideas of Ore and MacLane, J. Montes developed an algorithm to 
compute OM representations 

(1) t p = ('00, P i (01,p i ^l, P > Vh, P ) » ■ ■ ■ j (0r p ,p i , P ) VVp , P )! (0p ; ^ P > 0 P )) ) 

of each prime factor F p . An object t p as in m is a type', it contains several 
data structured into levels, encoding relevant arithmetic information about the 
polynomial F p and the prime ideal p. For the precise definition of a type, we refer 
to [8j. We now recall some of the properties of the invariants of a type. 

The number r p + 1 of levels is called the order of the type and r p is the Okutsu 
depth of Fp. The family of polynomials [0i lP ,..., (j>r v ,p\ is an Okutsu frame of F p . 
These are monic polynomials in 0[x] which are irreducible in O v [x]. 
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If we denote rm = deg^p, we have 

(2) mi | • • • | m rp | m rp+1 = n p , mi < • • • < m rp < m rp+1 = n p . 

The polynomial <f rp+ i tP := (f> p is an Okutsu approximation to F p ; it is a monic 
polynomial in 0\x\ of degree n p which is “sufficiently close” to F p for many pur¬ 
poses. More precisely, 

vOpW) > (Ap.p { 0 )) ■ 

m rp 

The data Ai, p ,A p are positive rational numbers called the slopes of the type. 
Typically, one denotes Ai,p = hj/ej, A rp+ i iP := A p = h Tp+ \/e rp+ \ the positive 
(coprime) numerator and denominator of the slope. One has e rp+ i = 1. 

The type determines a chain of finite extensions of the residue field F: 

F = F 0 cFiC...cF rp+ i=F p , 

where F p is isomorphic to the residue class field of the finite extension of K v de¬ 
termined by F p . The polynomial if o jP e F[yJ is one of the prime factors of the 
reduction of / modulo m. The residual polynomials f/'i.p e F,; [j/J are monic and irre¬ 
ducible; they determine the next extension: Fj+i = Fj[y]/(^j )P ), for all 0 < i < r p . 
The monic polynomial Vv p +i,p : = V’p 6 F p [y] has degree one. For i > 0, one has 
VAp ^ V- 

If we denote fi := degV’qp, we have 

(3) e(p/m) = ef-e rp , f(p/m) = f 0 ---f rp , m; = ei/i- • ■ 

We may consider each type t p as a path with root note V’o where each level is 
written along the edges: 

(,\ (<^l,p) Al lP ,l/>l I p) (^rp ,p ; A rp iP , t/Vp ,P ) (</>p, A p , ijjp) 

vV w o,p m . • tp 

Each node t of this path is identified with the type obtained by gathering all 
level data from the edges joining 4 with the root node. 

The last edge is dotted to emphasise that e rp +\f rp +i = 1, while for all other 
edges we have a ft > 1 as 0 and 0 show. 

Actually, the depth r p and the data (<fi, P , Aqp, ipi,p) of all levels i < r p are (up 
to certain equivalence relation) canonical data attached to F p [8] Sec. 3.3]. On the 
other hand, the last level (<j> p , X p ,ip p ) strongly depends on the choice of <j) p . 

For the proof of the next result, see Sec. 4]. 

Lemma 2.1. Let <fi e 0\x\ be any monic polynomial of degree n p that satisfies 
v ^ v (<j) p (6)). Then, for adequate choice of A e Z>o and a e F p , the object: 

tp = (V’o.p; (0i ,p 7 Ai,p) Vhjp); ■ ■ ■; ,p, A rp ,p, ipr t ,p)i (</>! A, y — a)) , 

is a type which constitutes an OM representation of F p too. 

The paths corresponding to the different OM representations computed by the 
Montes algorithm form a tree T of types. The leaves of the tree each represent one 
prime factor of / in O v [x\. The number of connected components of this tree (i.e. 
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the number of root nodes) is in one-to-one correspondence with the set of prime 
factors of / = / (mod m) in F[y]. For instance: 

■ • Tip 



In this example we have three prime ideals, which all share a common first level. 
Additionally, the prime ideals p and q share a common second level as well. Since 
the tree T is connected, / is a power of the prime polynomial ip$ in F[t/J. The 
polynomials F p , F q have Okutsu depth 2, while Fi has depth 3. 

Definition 2.2. For two prime ideals p, q e V, the index of coincidence 2 = ?'(p, q) 
is the first different level of their respective types. More precisely, 2 = 0 if 4 P and 
t q have different root nodes. If tfo.p — 00,q. then I is minimal such that 

(0£,P j Af jP , V’t.p) ^ (^,qi A^q, 0C,q)- 

One advantage of OM representations of prime ideals, is that they yield explicit 
formulas for the p-valuation of the ^-polynomials at each level of the type 4 q . The 
following proposition, extracted from |5] Thm. 3.1] and 0 Prop. 4.7], will be 
heavily used throughout the paper. It involves certain polynomials </>(p,q) e 0\x\ 
and certain hidden slopes A p , which are secondary data that have been conveniently 
stored along the running of the Montes algorithm [S Sec. 4], 


Proposition 2.3. Let p e V be a prime ideal of Ol■ Then for any 1 < i < r p + 1, 




Vi,p 

6l,p ’ ’ * l.p 


where Vy p = 0 and V+ i jP = e liP /y p (e ijP I0 p + ft i>p ) for all 1 < i r p . 

Let q e V be another prime ideal such that p A q and with index of coincidence 
I = i(p, q). For any 1 < i < r q + 1, 


Wp{<t>i,q(0)) 


' 0 , 

Vi + Xi 
ei ■ ■ ■ e-i- 1 ’ 

Vn + A p 

ei • ■ • e^_i ’ 

TUi, q Vi + Min {Ap, Aq} 
. mi e\ ■ ■ ■ ei-\ 


if 2 = 0, 
if i < 2, 

if i = 2 and <f>i A = 0(p, q), 
otherwise. 


In these formulas, we omit the subscript p. q when the invariants of the two types 
coincide. 


Okutsu bases of the integral closure of O v . Let L p be the completion of L 
with respect to the p-adic topology. We may consider a topological embedding 
L cz L p a K v , so that L p may be identified to a finite extension of K v of degree 
n p . We denote by O p the integral closure of O v in L p . 

Let r be the Okutsu depth of F p and suppose that 

It = (ifo, (01> ^1> 0l)i ■ ■ ■ i (0r + l) A r + 1 , 1pr+l)) ) 











6 


STAINSBY 


is the leaf corresponding to F p in the tree of an OM factorisation of /. 

The Okutsu frame [^> 1 ,..., cf r ] determines optimal polynomials <7o,p, • • •, g np -i,p 
in 0[x ] as follows. Each 0 ^ i < n p may be expressed in a unique way as: 

i = ao + a\mi + ■ ■ ■ + a r m r , 0 < ay < rrij+i/rrij = ejfj. 

Thus, the polynomials: 

fiTp : = x °' 0 He 0^i<n p , 

j =i 

are monic polynomials in 0[x] of degree degg^p = i. 

The following result is due to Okutsu El- 

Theorem 2.4. For all 0 ^ i < n p , the rational number w p ( gi, P (0 )) is maximal 
amongst all monic polynomials in O v [a:] of degree i. 

By Theorem ! 1. 41 we get a reduced triangular 0,,-basis of 0 P by taking the images 
under the embedding L a L p of the elements: 

(6) Oj := n~^ wp( ' 9i ’ p< ' e ^^gi }P (0), 0 ^i<n p . 

We call {oti)o^i <np the Okutsu basis of O p , or simply the Okutsu p-basis. 

If / is irreducible in O v [x ], then V = {p} and w = w p . Theorems 11.41 and 12.41 
show that in this case, © is a reduced triangular ^-integral basis of Ol- Thus, 
from now on, we may assume that ffV > 1. 

For further purposes, the family of numerators of an Okutsu p-basis is extended 
by adding an Okutsu approximation to F p , 

(7) A/p := {l = go,p> ■ ■ • j ffnp-i.p, 9np ,p := 4>p} ■ 

Definition 2.5. Let S c J> be a subset of prime ideals of Ol and consider the 
Okutsu set 



of all polynomials that are a product of exactly one extended Okutsu p -numerator 
for each p e S. 

By construction, all polynomials in Ok(5) are monic. Note that the Okutsu set 
depends on the choice of an Okutsu approximation <fi p as F p for each p e S'. 

Theorem 2.6. Let h e 0[x] be a monic polynomial of degree 0 < i < ns- For 
appropriate choices of the Okutsu approximations (f p , the set Ok(£>) contains a 
polynomial g of degree i such that 

w p (5(0)) 5 s w p ( h(8 )), V p e S'. 

This is one of the main results of the paper, whose proof will be postponed to 
Section [5j 

Up to finding the right Okutsu approximations, Theorems II.41 and 12.61 show that 
we may find a reduced triangular basis of Ol just by finding polynomials of degree 
0,1,..., n — 1, with maximal w-value in the finite set Ok("P) (c.f. Theorem 13.411 . 

Note that a brute force algorithm testing all possible factors < 7 i pjP of the Okutsu 
bases leading to polynomials of a fixed degree i would be exponential. 

A simple and very efficient algorithm, presented in Section [3j can be employed 
to choose an optimal combination of basis numerators for each degree i. 
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Okutsu bases of fractional ideals. Let I be a non-zero fractional ideal of O /,, 

I = p a > 
peV 

We consider a map giving a shifted valuation for a prime ideal p as a factor of 
the fractional ideal I: 

w pJ : L —>Qu {co}, 

a i—* w p j (a) = Wp (a) - a p /e(p/m). 

Also, for a subset S c p we define: 

w S,i ( ol ) = Min {wpj ( a ) : p e S} . 

Note that an element a e L belongs to I if and only if w-pj (a) ^ 0. 

Clearly these maps are consistent with the functions given in Definition 11.11 as 
Wp = w Pi o L and w = w V,O l - 

Definition 2.7. Let I be a fractional ideal of Ol- Let Sep be a subset of prime 
ideals of Ol and denote ns '■= 2 P es n p- An Okutsu S-basis of I is a triangular 
family 

go{9)/n^,--.,g na -mM v ^\, 

with numerators gi in Ok(£>) of degree i, such that Vi = wsj (9i(0)) > wsj ( h(6 )) 
for all monic polynomials h e 0„[;r] of degree i, and for all 0 < * < ns ■ 

A result analogous to Theorem 11.41 holds: an Okutsu iP-basis of I is a reduced 
triangular basis of / as an 0-module. Also, Theorem 12.61 holds if we replace w v 
with Wpj, because both functions differ only in a constant shift. Therefore, just as 
for the maximal order, it makes sense to compute an Okutsu S'-basis of I by looking 
for polynomials in Ok(S') with a maximal wgy-value amongst all polynomials in 
Ok(S') of a given degree. The MaxMin algorithm serves this purpose. 


3. MaxMin 


3.1. Formal extension of the Okutsu p-bases. The aim of the MaxMin algo¬ 
rithm is, given a set of prime ideals 5 c P and a fractional ideal I, to perform an 
efficient search for irsy-optimal polynomials in Ok(S'). 

To decide which numerators are chosen for each degree, we need only to know 
the values w q (gi p ^ p (9)) for all p, q £ S and 0 < i p =% n p . As presented in Section 
[21 these values are given by invariants present in an OM factorisation T of /. The 
exception is w p (4>p), which can be arbitrarily large, depending on the choice of <f>p 
the Okutsu approximation to F v . 

For this reason, we do not choose a concrete polynomial <p v beforehand, but 
rather run the algorithm as if Wp(<j>p(9)) (formally) takes the value oo. 


Definition 3.1. For all p £ S we define the following function on the Okutsu set: 


Wp : Ok(S) —* Q u {oo}, g 


(w p (g(6)), if(fp\g, 
[oo, if (f>p | g. 


This function does not depend on the choice of the Okutsu approximations cfp as 
by Lemma l2Jl and Proposition l2.3l the value of w q {<j>p{9)) for q # p only depends on 
p, q and not the choice of </> p . Thus, it makes sense to consider symbolic polynomials 
4>p of degree n p . 
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We consider a similar extension of the functions w p j, wsj to the Okutsu set: 
w pJ ( 9 ) : = w p (g) - a p /e(p/m), w S j ( g ) := Min {w p j (g) : p e S}. 

We have w$ i (g) < co for all g e Ok(S') with the exception of a single polynomial 
(f>s '■= ripeS^V Also, wsj (g) > wsj(g(0 )), and equality holds for adequate 
choices of all <f> p (depending on the given polynomial g e Ok(S'), g # </>$). 

The algorithm will provide a recipe to construct polynomials g t e Ok (S) of degree 
i with a maximal value of wsj ( gi ) among all polynomials of degree i in Ok(S'). 
The corresponding member of the triangular basis will be 

a.i = gi(9 )0 < i < ns- 

For a practical computation of on, we must apply the Single-Factor Lifting algo¬ 
rithm [9] to find concrete Okutsu approximations <p p , with a value w p (<p p (0)) large 
enough to guarantee that wsj ( gi ) = wsj ( gi(0 )) for all 0 < i < ns- 

3.2. The MaxMin algorithm. For each type t in the tree T, we denote by St £ S 
the subset of prime ideals p 6 S such that t is one of the nodes in the path joining 
the leaf t p with its root node. 

We fix an ordering S = {pi,..., p s } on the set S, with the property that for all 
types t in T, the subset S t c S is an interval of S. That is, there exist indices 
1 =% a± < &t ^ s such that, 

(8) S t = [a t , b t ] := {pj : a t < j < b t }. 

As the branches of T do not cross one-another, the reader will easily be convinced 
that it is always possible to consider such an ordering. 

We consider multi-indices i = ( i P ) P eS of degree degi := J] *p; leading to monic 
polynomials gi := Opffip.p e ^[xj in the Okutsu set Ok(S'), with degg; = degi. 

Definition 3.2. A multi-index i = (ip) pe s is said to be maximal if 

wsj (gi) ^ wsj (<?j), 

for all multi-indices j with degj = degi. 

In this case, we also say that gi is a maximal numerator. 

Notation. For 1 ^ j ^ s roe denote by uij the multi-index with coordinates = 0 
for all p # p j and i Pj =1. 


Algorithm 1 MaxMinfS 1 ] algorithm 

Input: A fractional ideal I of Ol and Okutsu numerators {gi tP : 0 < i < n p } of 
p-bases for each p e S. 

Output: A family io, ii,..., i„ s £ FF of multi-indices of degree 0, 1,..., n respec¬ 
tively. 

1: io <- (0, ...,0) 

2: for k = 0 —> ns — 1 do 

3: j <- Min {1 < i ^ s ■- w Pi j (g ik ) = w S j(g ik )} 

4: ife+i ife + u j 

5: end for 


The next result is the second fundamental result in the paper. Its proof will be 
given in Section [6] 
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Theorem 3.3. All output multi-indices o/MaxMin are maximal. 

This gives the name MaxMin for the algorithm, because it finds the maximal 
value amongst the minima of certain numerical data. This provides a computation 
of a reduced triangular basis as follows. 

Theorem 3.4. Let io, ii,..., i ns be an output o/MaxMin. Choose Okutsu approx¬ 
imations (j ) p of all p e V, such that 

wsj (Si fc ) = w S j (, g ik (60), 0 < k < n s - 

Then, gi 0 , ... ,gi n _ 1 are numerators of an Okutsu S-basis of I. 

Proof. Let h e (9„[a;] be a monic polynomial of degree 0 ^ k < ns- By Theorem 
12.61 there exists g £ Ok(S) (for adequate choices of all cf v for p e S) of degree k 
such that wsj ( g[9 )) ^ wsj (h(9)). 

On the other hand, regardless of the choices of the f> v , we have wsj(gi k ) ^ 
ws.i ( g ) by the maximality of ifc. Hence, 

WSJ ( gi k (60) = wsj (ffi J 5* ws.i (g) S* wsj (g(0)) ^ w S j ( h(0 )) . □ 

We will now present some remarks about the behaviour of the algorithm. We 
assume I = Ol for simplicity. 

3.2.1. Guaranteed termination. MaxMin always terminates after exactly ns itera¬ 
tions. 

Thanks to the convention Wp(<j > p ) = oo, the index j in step 3 indicates a prime 
p j such that for the multi-index ifc = (* p ) pe s, we will always have i Pj < n Pj . Thus, 
the next multi-index i/. +1 = (*p) p6 s constructed in step 4 has indices i' p < n p for 
all p. 

Furthermore, the first and last output multi-indices are io = (0,..., 0) and i„ s = 
(n P i,... ,n p J. As such, g io = 1 and g ins = Yl peS <l>p. 

3.2.2. Polynomial products are not computed. The algorithm does not compute the 
products gi k . It only computes the values w p (gi k ) for all p e S, which are determined 
by the 3-dimensional array of data w Vk (gj i .p i ) indexed by i, ji, and k in the ranges 
1 < * < s, 0 < ji =% n Pi , and 1 < k < s, respectively. 


3.2.3. MaxMin is not a universal maximiser. If the numbers w Pk (c/j ijPi ) are replaced 
by arbitrary, non-negative rational numbers Vk,j it i s Q>o and we take 


S 

i = 1 

with i = (ji)i^isgs a multi-index as above, the MaxMin routine may fail to compute 


Max {Min {vkj : 1 =% k < s} : deg i = d} , 


a maximal multi-index of degree d. 
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3.2.4. Initial conditions. Suppose i = (i p )p e s is a multi-index with degree deg i = d, 
such that wsj (< 7 i) is maximal amongst all multi-indices of degree d. Then, it may 
not be true that by increasing an adequate index by one, we get a multi-index j, of 
degree d + 1 , which renders a maximal value of ui(gj) amongst all multi-indices of 
degree k + 1 . 

For instance, let us consider the example presented in Section 13.31 The output 
index of degree 3 is 113 = (1,2,0), resulting in the polynomial g$ = </>i iP <(> 2 q with 
valuation vector w ( 33 ) = (18,12,12) for w p , and wi respectively. 

We could choose an alternative index j 3 = ( 1 , 1 , 1 ) which would give a polynomial 
<?3 = ( f > i,p4’i,q4 > i,i with the exact same valuations wig's) = (18,12,12). However, 
none of the indices (2,1,1), (1, 2,1), (1,1, 2) is maximal. For instance, J 4 = (1, 2,1) 
determines the polynomial g'± = <f>i t p(j)\„(j)i t \ with valuations w {g' 4 ) = (24,16,16). 
This is clearly not maximal as the polynomial constructed by MaxMin g 4 = (f>i t p(f> 2 ,q 
has valuations w ( g 4 ) = (18,22, 21 ). 

It is remarkable that the extremely simple strategy that MaxMin employs to 
choose successive maximal multi-indices is able to avoid these pathological cases. 

3.2.5. Ordering of input prime ideals. Theorem l3. 31 shows that MaxMin produces a 
sequence of maximal multi-indices regardless of the choice of ordering on S , as long 
as it satisfies ©. However, the numerators gi k produced from these multi-indices 
do depend on the choice of ordering. 

3.2.6. Complexity. To compute a triangular v-integral basis of L , a number of steps 
are required: 

(1) Use the Montes algorithm to produce an OM representation T of /. 

(2) Run MaxMinfP] to generate a family of maximal indices io, ..., i n -i- 

(3) Apply the Single Factor Lifting algorithm from [9] to get an adequate im¬ 
provement of the Okutsu approximation of each prime factor of /. 

(4) Compute the numerators of the Okutsu basis go, ■ ■ ■ ,gn -1 specified by the 
maximal indices. 

(5) Divide the Okutsu numerators by the appropriate power of 7 r to create an 
integral basis. 

The total complexity is equivalent to that of other OM based routines and is 
given in the following result. 

Theorem 3.5 ( [13t Thm. 6.1]). Suppose that F is a finite field with q elements 
and f is a separable polynomial. Take 5 := u (disc(/)). The total cost of the 
computation of a v-integral basis of L by the application of the Montes and the 
MaxMin algorithms is 

O ( n 2+e S 1+e + n 1+e S\og(q) + n 1+e S 2+e ) , 

operations in F. If we assume q small, this will give us a refined estimation of 
O (n 2+£ <5 1+£ + n 1+£ <5 2+£ ) bit operations. 

3.3. MaxMin example. We will now present a small example for I = Ol and 
S = V = { p, q, l}, where T is connected. The tree is shown in Figure [U where we 
indicate only the data (c t> , A) for each edge. 
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Figure 1 . Example connected tree T of types. 


Since all slopes have integer values, all denominators ei are equal to one. We 
assume that /o = mi = deg ipo = 1 and: 

p : ei = l,/i = 4, hi = 6 ; 
q : ei = 1, /i = 3, fti = 4; e 2 = 1, = 2, /12 = 6 ; 

1: ei = l,/i = 3,/n =4. 


Note that n p = 4, = 6, and rq = 3, so that n = 13. 

The data corresponding to the edges leading to a leaf are not specified as we do 
not need them to run MaxMin. 

Suppose moreover that 

0(p, q) = Hp, 0 = H,p = H.q = </>i, 1 = 0, 0(q- 0 = fa, q = 0', 

and the hidden slopes are: 

A2=A; = 6, AP=Af=4, A; =6, A? = 5. 

The numerators of the extended Okutsu bases of each prime ideal will be, 

: 1) 01,p) 01, p, 01, p, 0pi 

•A/q : 1, 01, q, 01, q, fa, q, 02,q01,q; 02,q0l,qi 0q > 

M : 1, 01,i, 01,[, 01- 

Using the explicit formulas of Proposition 12.31 we may compute the valuations 
of each of the ^-polynomials. We write them as a tuple w = (w p ,Wq,Wi). 

w(0i, P ) = (6,4,4), w(fa) = (oo,16,16), 

wfa hq ) = (6,4,4), w (02,q) = (12,18,17) w fa,) = (24, 00 , 34), 

w(0i,[) = (6,4,4), w(0i) = (12,17,oo). 


We can now step through the results of running MaxMin. The “minimal” valua¬ 
tion is underlined at each step. This indicates the index which will be incremented 
in the following step. 


i 

9% 

w(9i) 

w{9i) 

0 

1-1-1 

(0,0,0) 

0 

1 

01, p • 1 • 1 

(6,4,4) 

4 

2 

01,p * 01, q ’ 1 

(12,8,8) 

8 

3 

01, P • 01,q • 1 

(18,12,12) 

12 

4 

01,p ' 02,q ' 1 

(18,22,21) 

18 

5 

01,p ‘ 02,q • 1 

(24, 26, 25) 

24 

6 

01,p ’ 02,q ' 1 

(30,30,29) 

29 
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7 

01,p ’ 02,q ’ 01,1 

(36,34,33) 

8 

01,p ' 02,q ' 01,[ 

(42,38,37) 

9 

01,p ' 02,q ' 01 

(42,47,oo) 

10 

0p ’ 02,q ’ 0[ 

(oo, 51, oo) 

11 

0p ' 02,q01,q ' 01 

(oo, 55, oo) 

12 

0p • 02,q0i,q ■ 01 

(oo, 59, oo) 

13 

0p ' 0q ' 01 

(oo, oo, go) 


The final element g 13 is not included in the v-integral basis. 


4. Computational examples 

In this section, we will present a number of example computations using an 
implementation of the MaxMin algorithm for the computer algebra system Magma 
[2]. We compare MaxMin’s execution time for computing u-integral bases with that 
of the method of the quotients, another OM-based algorithm, as well as the internal 
routines found in Magma. 

All executions were performed on GNU/Linux running on 8 -core 3.0GHz nodes 
with 32GB main memory. Each execution ran in a single core, using Magma 2.18-5. 

Examples will be given for number fields L = Q[x]/(f) for polynomials / e 
Z[x\ and function fields L = F q (t)[x]/(f) over a finite field F g for polynomials 
/ e F g [f] \x\. The example defining polynomials are taken from [9]. 

The first example is comprised of the H-class of polynomials, 

B p ,k{x) = ( x 2 - 2x + 4) 3 + p k , 

of degree 6 . We take /( x) = i?i 3 ,fc(:r) e Z[x\ in the number field case and f(x) = 
Bt3+2,k(x) e F 7 [t, x] in the function field case, with k < 5000. The execution times 
for computing a Hermitian u-integral basis of L are shown in Figure [2j 




k k 

Bi 3 ,k(x) e 7L\x\. 

Bt3+2,k (x) e F 7 [t,x\. 

Figure 2. Running time for maximal order Hermitian p -basis 
computation defined by polynomials B p k(x) with k < 5000. 


Due to the low degree of the field extensions, computing the HNF of a basis is 
negligible compared to computing the basis itself. It can be seen that both OM- 
based routines have a similar performance and are roughly 100 times faster than 
the internal Magma routines in the number field case for k = 5000 and 1000 times 
faster in the function field case. 
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In order to demonstrate the performance of MaxMin on larger polynomials, we 
will consider the H-class of polynomials 

A£ n>k (x) = ( x n + 2p k ){(x + 2) n + 2 p k ) ■■■((x + 2 m- 2) n + 2p k ) + 2p nmk . 

These polynomials have degree nm. In the number field case, we take f(x) = 
Af 0 i 29 ( x ) e %[x] with n < 100. Figure [3] shows the times for the OM-based 
methods and the total times when we include the time to compute the Hermitian 
basis. 




Defining polynomial degree Defining polynomial degree 

Without HNF. With HNF. 

Figure 3. Running time for maximal order p-basis computation 
defined by polynomials Af 01n2 g(x) with deg(H| 01 29 ) ^ 400. 


From this example, we can see that MaxMin is somewhat faster than the Method 
of the Quotients, however when the time to compute the HNF of the resulting basis 
is included, we see the advantage of the triangular basis computed by MaxMin. In 
this example, Magma took 257 seconds to compute the basis for deg / = 40 (and 
does not appear in the figure), and was unable to complete the computation for 
deg / = 80 due to main memory limitations. 

We consider a slightly smaller example in the function field case, with f(x) = 
4L 2 „6 e F 7 [t,ir], once again for n < 100. The resulting execution times are 
presented in Figure |U 


90 
80 
«T 70 
§ 60 
~ 50 
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10 
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Figure 4. Running time for maximal order p-basis computation 
defined by polynomials Al 01 n 29 {x) with deg(H| 01 29 ) =% 400. 
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Here, MaxMin computes a basis considerably faster than the Method of the 
Quotients, however the improvement in computing the HNF from the triangular 
basis is not as pronounced as in the number field case. The second data point 
for the Magma routine is at deg / = 60, which took 3304 seconds to compute a 
Hermitian basis. 

Finally, we will consider some very large examples. We take the recursively 
defined AC-class polynomials defined in this case as EC p j(x) = E p j(x)-C Py 28 +P 900 , 
where the A-class and C-class polynomials are given in [5]. Firstly, we consider 
f(x) = ECioi y s(x) e Z[x], a degree 900 polynomial which splits into six degree 6 
factors and one degree 864 factor over the p-adics. Secondly, we consider f(x) = 
EC t 2 +4A (x) e F 7 [t, re], a degree 72 polynomial that splits into four degree 9 factors 
and one degree 36 factor over the p(£)-adics. 

Table 1. Running time (in seconds) for maximal order p-basis 
computation defined by polynomials EC Py j. 


ECioi,s(x) e Z[x], EC t 2 +iA (x) e F 7 [i,a;]. 


Algorithm 

Basis 

HNF basis 

Algorithm 

Basis 

HNF basis 

MaxMin 

9.9 

112.6 

MaxMin 

13.3 

21.5 

Quotients 

21.1 

429.3 

Quotients 

89.5 

8353.8 


These final examples show advantages of using MaxMin to compute a local basis, 
whether or not a subsequent step to pass the basis to HNF is required. However, 
they also illustrate the advantages of using a triangular basis instead of a basis in 
HNF where it is possible. 

5. Optimal polynomials in the Okutsu set 

In this section, we will present the additional details necessary to prove Theorem 
12.61 This will be broken into two parts, Proposition l5.16l and Proposition l5.231 which 
together prove the theorem. 

5.1. Operators associated with a type. Consider a type t of order r over (A, v): 
t = (i/> 0 ; (</>i, ...; (0 r , A r ,i> r )) . 

The type 4 carries three kinds of operators. There are normalised valuations: 

Vi : A' t ,[a;] —> Z u { go } , 0 =% i < r. 

The last valuation v r is also denoted v t . Also, we have Newton polygon operators: 
Nj = N w _ 1;0i : K v \x\ —* 2 r2 , 1 < r, 

where 2 R is the set of subsets of R 2 . The image of 0 is the empty set. For every 
non-zero g e I\\x\ we consider the canonical 0 j-expansion g = So < s a s ( t > ii where 
a s e K[x\ have degree less than deg0j. Then N i(g) is the lower convex hull of the 
set of points (s,Vj_i (a s 0 f)). 

Finally, we have residual polynomial operators: 

R-i = '■ K v[x\ —* Fj[y], 0 ^ i < r, 

which are multiplicative: R i(gh) = R^g) Rj:(Ii) for all g,h e A'[x]. 
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The valuation vq is defined as 


vo I a s X s j = Min {v (a s ) : 0 < s} . 

\0«: s / 

The residual polynomial operator Ro is defined as: 

Rote) = 7r - " o(s tete) (mod m[y]). 

For 1 < * < r the operators Vi, Ri are defined in a certain recurrent way (SJ Sec. 

3-1]. 

Definition 5.1. We say that t is optimal if mi < m 2 < ■ • ■ < m r . 


Definition 5.2. Let t be a type of order r 3 s 0 and g e K v \x\. We define ord t (g) 
as the non-negative integer ord^, r (R r (g)) where ip r and R r (g) are considered as 
polynomials in F r [y]. If or d t (g) > 0 we say that 4 divides g and we write t \ g. 

Definition 5.3. Let t be a type of order r and g e 0\x\. We say that g is a 
representative of t if g is a monic polynomial of degree m r+ \ := e r f r m r such that 
R r te) = ipr- We denote by Rep(t) the set of representatives oft. 

If <p £ 0\x\ is a representative of t, then by definition t = (t; (<p, A, ip)) is a type 
of order r + 1 for any choice of a positive rational number A and a monic irreducible 
polynomial ip e F r +i[y], ip # y- 

Definition 5.4. A prime polynomial is a monic irreducible polynomial in <!?„[a:]. 

Since R r is multiplicative and tp r is irreducible, the representatives of t are prime 
polynomials. 

The next result collects the essential properties of the polynomials which are 
divisible by a type. It is taken from (SI Lem. 1.4, Thms. 3.1, 3.7], 

Theorem 5.5. Let t be a type of order r ^ 0 and let <p e 0\x\ be a representative 
oft. Let F e 0„[x] be a prime polynomial and choose 6 e K v a root of F. 

If <p A F and t | F, then 

(1) N Vr ^{F) is one-sided of slope —A, for a certain positive rational number A 
such that v (<p{9)) = (V r+ i + ^)/( e l ’ ' ‘ e r)- 

(2) deg(F) = deg(^), where I = I(N Vrt ^(F)) is the length of the Newton 
polygon; that is the abscissa of the right end point. 

(3) deg(R„ rj 0 i A(I ;1 )) = and R. Vr ,<f>,\(F) = ip a , where e\ is the least positive 
denominator of X and ip e F r+ i[y] is a monic irreducible polynomial ip A y. 

(4) The type t = (t; (ip,\,ip)) divides F. 

For 0 < i =% r, let t, = Truncj(t) be the type of order i obtained by dropping all 
levels (<pj,\j,ipj) of order j > i from t. 

By the definition of a type, each <pi is a representative of tj_i, for 1 ^ ^ r. 

The next result follows easily from Theorem 15.51 and J5J Lem. 2.4]. 

Theorem 5.6. With the above notation, tj | F for all 0 < i ^ r, N i(F) is one 
sided of slope —A and R,(F) is a power ofipi. 

Finally, the following result follows from Theorem 15.51 and [5] Prop. 3.5]. 
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Theorem 5.7. With the above notation, let g e 0„[x] be another prime polynomial 
such that t | g, t \ g. Then, 

, taw de §(5) V r +1 + Min A'} 

•’ <9<< ' )) “*iW-Spry;-• 

where —A' is the slope ofN Vrt< f,(g). 

5.2. Non-optimised tree of types. Let us briefly describe how the Montes algo¬ 
rithm constructs the tree T of OM representations of the prime factors of the input 
polynomial / e 0\x\. Initially, / is factorised in F ['</]. For each monic irreducible 
factor ip of /, a triplet (t, </>, w) is considered, where t = (<p) is the type of order 0 
determined by ip, (f> is a representative of t (that is, a monic lift of ip to 0 \x ]), and 
w = ord t (/) = ord v (/). All these triplets oj) are stored in a stack. 

Along the execution of the algorithm the stack always contains triplets (t, </>, w), 
where t | /, <fi is a representative of t and u> = ord* (/). The main loop of the 
algorithm takes such a triplet and attaches to the type t one or more branches 

,if, '■= (t; (</>, A ,ip)) of t such that \ f and the pairs (A ,ip) are considered as 
follows, 

• —A runs on the slopes of N"^/) := N„ t> 0 (/) n ([0,w] x R). 

• runs on the prime factors of R« t , 0 ,A(/)- 

Let cj) A,-i/> e @[ x ] be a representative of tA,^ and take wa,^ = ord^,(R l , t> 0 > A(/)) = 
ord tA ^(/). If this positive integer is equal to one, then divides only one of the 
prime factors F v of / in O v \x\. In this case, we add a final level to t to construct 
the leaf t p associated with this prime factor. On the other hand, if ui\^ > 1, then 
the triplet </>a,Vo ^a,^) is pushed back onto the stack to bare further branching 

in future iterations of the main loop. 

After a finite number of iterations of this process, the algorithm outputs a list 
til,..., tiv of types parametrising the prime factors of / in O v \x\. Let us denote by 
T nop the tree of types obtained by this procedure. Note that for every node t e T nop 
which is not a leaf, the edges with left end point t have the same ^-polynomial; a 
tree of types with this property is said to be coherent. 



Figure 5. A segment of the non-optimised tree T nop . 

This describes a kind of “non-optimised” Montes algorithm, yielding a “non- 
optimised” tree of types. The types it a,v> ma y n °t be optimal. In fact, if A e Z and 
deg ip = 1, we have 


deg (j> a,v> = e A • deg • deg <\) = deg <£, 
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where e\ is the positive denominator of A. We must avoid this situation, because 
the numerical data attached to the types will not be intrinsic data of the prime 
factors of /. 

For this reason, the Montes algorithm includes a “refinement procedure” which 
ensures that it only stores optimal types. However, a price must be paid; the output 
tree of OM representations is no longer coherent. 

The optimised tree of OM representations (the real output of the Montes algo¬ 
rithm) may be derived from the non-optimised tree by an iterative application of 
the following transformation. Any path, 


(9) 


(</> i,Ai,^i) 

H •.• 


(<^71) An, lf n ) ,/ 

•-• 11 


in which all edges except for the final one are bad edges satisfying A; G Z, deg ift = 1 
for i < i < n, collapses into 

(10) t . * t „ with A* = Ar + • • ■ + A„. 


The types t' and t" are “equivalent”, and this means that P(t) = P(ti'), where P(t) 
is the set of prime polynomials g G O v [x] which are divisible by t [121 Thm. 3.7]. 

The existence of the non-optimised tree is useful in many situations. Let us see 
an example. 

Lemma 5.8. Let t,t' e X nop . If t is a truncation oft', then P(t) z> P(t'). If 
neither of these types is a truncation of the other, then P(t) n P(t') = 0. 

Proof. The first statement is an immediate consequence of Theorem 15.61 

The second statement is obvious if t and t' have different root nodes, because 
for all F G P(t), the reduction F modulo m is a power of the rnonic irreducible 
polynomial if o corresponding to the root node of t. 

Suppose that t, t' have the same root node and let t 0 be the greatest common 
node in the paths joining t, t' with their root node. By the first statement we may 
assume that t and t' are branches of to, in other words, that to is the previous 
node of both t an t'. By the coherence of T no P we have 

t = (t 0 ; (</>, A, if)), t' = (t 0 ; (</>, A', if')), 

where either A ¥= A' or A = A' , if / if' ■ 

Let r be the order of to and v r its attached valuation. Now, for any F G P(t), 
F' e P(t'), Theorem 15.61 shows that N Vr ^{F) and N Vr ,^(F') are one-sided of slopes 
—A and —A' respectively. Hence, A ^ A' implies F # F'. On the other hand, if 
A = A' then Rt o , 0 ,A (F) = if a and Rt o , 0 ,A(-F , ) = ( if') a and this implies F / F', 
because if if'. □ 

This result may be false for arbitrary incoherent trees. However, Lemma |5.81 is 
valid for the optimised tree ‘X of OM representations of the prime factors of /. 

Proposition 5.9. Let t,t'el he two nodes such that neither of them is a trunca¬ 
tion of the other. Then P(t) n P(t / ) = 0. In particular, Vt n TV = 0• 

Proof. Clearly, the nodes t,t' are equivalent to two nodes of the non-optimised 
tree, neither of them being a truncation of the other. Thus, the statement is an 
immediate consequence of Lemma 15.81 

The final statement is a consequence of Vt = {p e V '■ F v G P(t)}. □ 
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Consider the chain of refinements that take place between © and m- During 
each refinement that provokes branching of a type, the intermediate <j> and A values 
are stored. 

Let tp,tq G T be two leaves attached to prime ideals p and q with index of 
coincidence i (p, q) = £. Then suppose that at level £, each type has a list of stored 
refinements, 

Refute,) = ■ 

This allows us to extend the index of coincidence to a more precise indicator. 


Definition 5.10. The minor index of coincidence z(p, q) for two leaves t p ,t q G T, 
is the least index t', such that for the refinement lists given in ED, 





We also define the extended index of coincidence of two types as, 


/(p,q) := [z(p, q), z(p, q)]. 


These extended indices of coincidence are ordered lexicographically. 


Definition 5.11. Let t p , t q G T be two leaves with index of coincidence z(p, q) = l 
and let the list of refinements of each type at level t be as in ED- 

(1) The greatest common ^-polynomial of the prime ideals p, q is q) = 
<f> J) = <fi,f, with j maximal. 

(2) The (non-optimised) hidden slopes of the prime ideals p,q are (A p ) nop = 
A*Jj and (A q ) nop = AH, for this maximal value of j. 

(3) The (opimised) hidden slopes of the prime ideals p, q are A p = A^ + - • 
and A q = A*^ + • • • + A*^, for the same maximal value of j. 

Remark 5.12. (1) By (flOl) A^ p = particular A^ jP ^ A p , for all qeT 

with z(p, q) = l. 

(2) Proposition 12.31 is easily deduced from Theorem 15.71 Since this theorem is 
valid for arbitrary types, it is clear that the formulas in Proposition l2.3l are valid for 
the (^-polynomials of the non-optimised tree just by replacing the optimised hidden 
slopes with the non-optimised ones. 


5.3. Optimal polynomials as products of (^-polynomials. Let S <= V be a 

subset of prime ideals. Let T s be the tree of OM representations of the prime 
ideals p g S computed by the Montes algorithm. We keep the content of Section [2] 
concerning the data attached to the different types t p for p G S. We recall that the 
polynomials G 0\x\ are concrete choices of Okutsu approximations to the prime 
factors F p of /. 

The (^-polynomials for all the prime ideals generate a semigroup. 

Definition 5.13. Let be a set of prime ideals. We denote by $(5) c: 0[ x] 

the multiplicative semigroup generated by 

{(fi.p : p g5, 0 < i ^ r p | u (J App(p), 

p eS 
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where App(p) = {<j> £ 0[x] monic of degree n p such that w p ^ w p (</> p (#))}. 

We use $(p) to denote $({p}). 


We are interested in showing that we can restrict our search for polynomials of 
a given degree d with maximal ws,/-value to those in the semigroup < h(S'). 


Definition 5.14. Let g £ 0[a;]. The degree adjusted w p -valuation of the element 
g(9) e Ol is defined as 


™P : = 


Wp (g(°)) 

deg g 


Lemma 5.15. Let t be a node in the non-optimised tree T nop and let g,h e O 0 [x] 
be two prime polynomials divisible by t. Then, for any prime ideal p e V\Pt we 
have w p {g(9)) = w p ( h(9 )). 


Proof. If t and t p have different root nodes, we have w p ( g(9 )) = 0 = w p ( h{6 )), 
because F p is a power of tf o. p and <?, h are powers of the root node of t. 

If t and t p have the same root node, let to be the greatest common node in the 
paths of T nop joining t p and t with the root node. Since p f V t , the node to cannot 
be equal to t. Since t p is a leaf of the tree, to cannot be equal to t p either. The 
structure of the non-optimised tree is shown in Figure [6] 



Figure 6. The node to is the greatest common node of t and t p . 


Let t',t" be the nodes following to in each of the two paths. Since the non- 
optimised tree is coherent, we have 


t' = (t 0 ; (0, A', V’ , )) 5 t" = 


with a common choice for the representative </> of to- 

By Theorem l5.61 N v t0 ,<f>(g) and N Vt0 , 4 >{h) are one-sided of slope —A". Proposition 
15.91 shows that t' Jf g, t' Jf h. By Theorem 15.71 we have 


wp (g(9)) = 
where r is the order of to- 


V r+1 + Min {A', A"} 


deg cj) 


ei 


= w p (h(0 )), 


□ 


The next result is the main aim of this section. 


Proposition 5.16. Let S c J> be a set of prime ideals. For any h e 0^[a;] monic 
of degree 0 < d < n, there exists (f> e <&(fP) also of degree d such that, 

(12) w p (0(6*)) 3s w p (h(9 )), V p £ S'. 

Proof. The proof will proceed by induction on the degree d of the polynomial. We 
will work in steps, in each one reducing the space in which we need to consider h. 
If d = 0, then <f> = h = 1 £ ^(S). 
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Claim. It is sufficient to check (11211 for h a prime polynomial. 

Let h = ftift- 2 , with hi,h -2 e 0„[a;] monic of degree di,d ,2 > 0 respectively. 

By the induction hypothesis, there exist 4> t e $(S I ) of degree di such that, 

w p > w p (hi(9)), V p e S, 

for i = 1, 2. Then, <j> = (f>i4>2 £ 'I* (•S') clearly satisfies (fT^l) . This proves the claim. 

Now, assume that ft, is a prime polynomial. If gcd(/, ft) = 1, then w p ( h(6 )) = 0 
for all p £ S. Thus, (fl2l) is obviously satisfied. 

Therefore, we can assume that ft = i/jq, kN, for i/jq e F[y] a prime factor of /. 
By hypothesis, the root node ipo (thought of as a type of order zero) divides ft. 
Let t be the highest order node in the non-optimised tree Tg° p such that t | ft, and 
let i be the order of t. We distinguish two cases according to t being a leaf or not. 

Case 1. t is a leaf. 

In this case, S t = {po} contains only one prime ideal. The ^-polynomial in the 
last level of t is an Okutsu approximation <fi Po to F Vo . Since t | ft, Theorem 15.51 
shows that deg h = I ■ deg (f> Po = I • n Po for some positive integer £, and z> po (h(9)) > 

v Po OM 0 ))- 

Let us consider (j) 0 e App(p 0 ) and close enough to F Po so that 

w Po (MO)) > u> Po (h(6))/£, 

and take <t> = (j> Je $(iS). By construction, ui Po (<t>(6)) > w Pa (h(6)). On the other 
hand, for any p e S, p A po, we clearly have w p (</>(#)) = w p ((/>o{0)) and Lemma 
15.151 shows that w p (4>o(9)) = Wp (h(0)). Since deg(<^) = deg(ft) we deduce that 
w p (cj>(0)) = w p (h(6)). This proves (fT^l) . 

Case 2. t is not a leaf. 

For a certain choice (f> t of a representative of t, the node t has several branches 
in the non-optimised tree, of the form 

("t5 (^t) ^5 VO)- 

By the maximality of t, we have \ ft for all these branch nodes. Let A max be 
the greatest slope (in absolute size) of these branches and let t max be any branch 
node of t with slope A max . 

Since t | ft, Theorem 15.51 shows that N t , t ,<p t (ft) is one-sided of slope —A h and 
deg(ft) = ^deg((/> t ), for certain positive A^ G Q, £ e Z. 

If for some po £ St we take </>y Po satisfying: 

(13) deg(</>j, Po ) = deg(4> t ), * I <t>j, Po> 

then Lemma [5.151 shows that w p (h(6)) = iu p (<j>j,p 0 (0)) for all p f S t . As in Case 
1, for (j) = e ®(S) this implies w p ( h(9 )) = w p (<f>(0)) for all p £ S t . Therefore, 
we need only to find some (f>j, Po satisfying (fl3l) and w p (<j>j t Po (0)) ^ £~ 1 w p (h(0)) for 
all p G S t . Then we shall have da. 

Let t' be any node of the optimised tree which has been derived from t max by a 
series of refinement steps as indicated in © and (fTU|) . 

Let (<(>o, Ao, ipo) be the last level of t' in the non-optimised tree. As explained in 
Section 15221 the last level of t 1 as a type from the optimised tree will be (<j> o, \q,iJjo), 
where Aq is the sum of all the slopes of all bad levels between t' and its previous 
node in the optimised tree. 
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^max 



(00 Aoj V'o) 

•-• t' 


Figure 7. The node t' corresponds to a node of the optimised tree. 


Thus, 0o = 0j,p o for all po e -St', where j is the order of t' as a type of the opti¬ 
mised tree. Clearly, 0j, Po satisfies (THU) : let us compare w p (0j, Po (0)) and w p ( h{6 )) 
for p e S t - Take p e S t and let tx,^ be the unique branch of t such that \ F p . 
By Theorem 15.71 

d —i futaw ^+i + Min {Aft,, A} V i+ i + A 

i Wp (h( 9 )) = - -. 

ei • • • e\ ■ ■ ■ ei 


Thus, we need only to show that 


Wp (00(0)) S* 


Vi+l + A 

€i • ' ’ €{ 


Suppose that t max gives rise to a node of the optimised tree. In this case, we 
have 4 max = t' and 0 O = 0 t . By Theorem[531 w p (0o(0)) = (Vi +1 + A)/(ei • • • ef). 
From now on we assume that t max / tib 

If It^ ^max, then t | 0o, t\^ Jf 0o by Proposition 15.91 hence Theorem 15.71 
shows that 


. . .. _ V)+i +Min{A , Amax } 

w P (0o (0)) = - 

ei ■ ■ ■ e, 


f0+i + A 

ei • • • 


Finally, suppose that tx,ip = t max . The order of t' as a type of the non-optimised 
tree is ^ i + 2. By Proposition 12.31 applied to the non-optimised tree (see Remark 
15.121) we have 


w p (0o(0)) > 


Vi 


+2,p 


ei+i,p/j+i,p( e i+i,pTj+i + hi + i, p ) _ I0+i + A 


Cl • • • GiCi- 1-1,p Cl p Ci • • • 6^ 

because ei+i, p = /j+i, p = 1 and A = /q+i, p . Thus, in all cases we obtain the desired 
inequality. □ 


5.4. Optimal polynomials as products of numerators of Okutsu bases. By 

Proposition 12.31 for any p e S we have 


(14) w p (0q p (0)) 


1 Vi + Xi 

777-2 Cl ' • * C 2 — 1 


1 V i+ i 

777-2+1 Cl ’ C 2 


< Wp (0z + !,p (0)) , 


1 Sj * ^ 


P- 


Let us analyse how closely we can replicate this inequality (fill) for cross valua¬ 
tions, that is to say when the 0-polynomial belongs to a different prime to that of 
the valuation. The next results follow closely form Proposition 12.31 too. 


Lemma 5.17. Let p, q e V be two different prime ideals with index of coincidence 
l = i{ p, q). Then: 

(1) Wp ( 02 , q (0)) < Wp (02+1, q (0)), 1 < i < i, 

(2) Wp (02, q (0)) = Wp (02+!, q (0)), t < i < r p . 
















22 


STAINSBY 


It is easy to find examples where 

(15) w p (<&, q (0)) > w p {<j> t+ i, q (6»)). 

This pathology occurs when (f>e,q = <fi(p, q) and Ap is much larger than Aq (see 
Proposition 12.3|) . Hence, it is also easy to find specific conditions that avoid (fl5l) . 

Lemma 5.18. Let p,q e V be two different prime ideals with £ = z(p,q) > 0, 
chosen so that Aq 0 Ap. Then, 

(16) Wp (<^,q(0)) = Wp (<f>e+i, q {0)). 

In particular, every numerator of degree i of the Okutsu q -basis has maximal 
p -valuation amongst all polynomials <j> e $(q) of degree i. 

Proof. By the hypothesis, Min {Ap, Aq } = Ap and Proposition ^. dl gives (1T51) . There¬ 
fore, Lemma 1 5.1 71 and (1161) show that the p-valuations of the polynomials bi.q in¬ 
crease with their degree up to index £ and then remain equal. As such, we will 
always have a maximal valuation by taking higher degree ^-polynomials, rather 
than products of smaller degree ones. □ 

Using the extended index of coincidence presented in Definition 15. 101 the follow¬ 
ing Lemma gives us a link between the relative similarity of prime ideals and their 
respective cross-valuations in certain cases. 


Lemma 5.19. For a prime ideal q e V, let p, [ e P\ {q} be two prime ideals such 
that either l = p or they satisfy: 

(1) 1 (1, q) ^ J(p,q), and 

(2) A?>A< if 1(1, q) = J(p, q). 

Then, for £ = i(\, q) and me = = me,i, we have 

wp (^™' l ’ l ' /me (9)^ ^ Wp (0i, q (0)), £ i r q + 1. 


Proof. We may consider four cases depending on the relationship between the three 
prime ideals p, q, and 1 . 


Case 1. i(\, q) > z(p,q). In this case, k := z(p, l) = z(p,q) < £ and Aq 
Ap = Ap. Therefore, for £ < i < r q + 1, Proposition 12.31 shows that 


(17) 


Wp 


q t m t 


wii.q meVk + Min {Ap, Af } 
m k Cl ■' ■ C/e—1 


mj.q 14 + Min (Ap, Aq} 

TYl k 6i Cfc—1 


Wp (<t>i, q(6>)) . 


Af 


[ ! 


Case 2. z(l, q) = z(p, q) and either l = p or z(p, () > £. We have 


m iiq /mi 

‘e,i 


rrij q Vi + Xe,p 
me e\---ee-\ 


by Proposition [331 On the other hand, 
' V e + X q p 

(18) Wp (4>i,q(9)) = < 


e\ • • ■ ee -1 

TO i,q Vi + Min (Ap, Aq} 


if i = £ and 


otherwise. 


me 


ei ■ ■ ■ ee-i 


^(p,q), 
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By the first item of Remark 15.121 we have A p < A^ p , so that w p ((/>,, q (0)) < 
w v (^Ti q/,mf (^)) i n both cases. 

Case 3. «([, q) = *(p,q) = f(p, l), and 1(1, q) > /(p,q). In the non-optimised tree, 
we find the situation described in Figure[8](a), where we have written the optimised 
hidden slopes instead of the non-optimised ones. 

We necessarily have </>^ q ¥= </>(p,q) = </>(p, 0 # (j>t,u and Ap = A p , A q = Af. 
Therefore, the equations of (11711 (where k = £ now) are again a consequence of 
Proposition 12.31 for all i < i < r q + 1. 



Figure 8 . Relative positions of t[, t q , and t p in the non-optimised 
tree when i(l, q) = i( p, q) = i(p, l). 


Case 4. p ¥= 1, z(l, q) = i(p,q) = i(p,l)j and 1(1, q) = /(p,q). In the non-optimised 
tree, we find the situation described in Figure [8] (b). We have </>(p,q) = </>(!, q) = 
</>(p, 1). By our assumptions, Ap < Af and Proposition 12.31 shows that: 


w p 



q Vi + A|, 

me ei ■ ■ • e^_i ’ 


whereas w p (</>i, q (0)) is given by (fl8l) . Since Ap = Aj,, we get w p (^i>; :q ( 9)) < 
W P ’ q ^ me (0)j as desired. □ 


Deflnition 5.20. Let g = ]^[ p G The disorder of g is calculated as 

D (s) = 2 P eS max {deg(^j p ) - n p ,0}. 

Definition 5.21. A polynomial p p = ^“p^ip • • • </Vp( p <?! > p r, ’ +1 e $(p) sa ^d to he 
canonical if 0 < dj < ej )P /j iP for all 0 ^ i =% r p . 

The canonical polynomials G d^p) of degree deg(y> p ) ^ n v coincide with the 
numerators of the Okutsu p -basis, and hence belong to Ok(p) too. 


Remark 5.22. For a given set of prime ideals S ^ V, all of them with the same 
root node, it is always possible to choose a prime ideal po e S such that Ap 0 < A p ° 
for all p G S. 

To do so, begin at the root of the tree of types representing S and move up 
through the levels. When branching is encountered, take the branch that corre¬ 
sponds to the slope of least absolute value at that level. We continue in this way 
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until we reach a leaf node, which will correspond to a prime po with the desired 
properties. 

Proposition 5.23. Let S ^ V be a set of prime ideals and consider <j> G <I>(S') 
monic of degree 0 =% d ^ ns- For appropriate choices of the Okutsu approximations 
<pp , the set Ok(5') contains a polynomial g of degree d such that, 

w P (9(0)) > wp 0(0)) i V p e S. 

Proof. Consider <j> = Opes Pp the separation of the polynomial 0 into its p-parts 
tpp for each peS. Then, let So = {p e S : <p p =£ 4>p}- 

We will follow an iterative sequence of three steps to find a polynomial g 6 Ok(S) 
that meets the requirements of the proposition. Throughout this process we will 
be modifying g, which is initially set to 0 , via its individual p-parts: 

(1) For all p G So, make <p p canonical. 

(2) If D(g ) = 0, then for all p G S\S 0 we take tp p = 0 p to be an Okutsu 
approximation to F v with w p (0 P (0)) ^ w v (0(0)), and finish the iterative 
process. 

(3) Fix some q G So with deg(</ 3 q ) > n q and then select [ e So, the “closest” 
prime ideal to q. Transfer all ^-polynomials in <p q , except for a single 0 q , 
to ipi. Remove q from So- 

Below, we will show that at each step, the p-valuation of g does not decrease 
for all p G Sq and that the process will terminate in a polynomial g belonging to 
Ok(S) after a finite number of iterations. The inequality w p (g(9)) > w p (0(0)) for 
the primes p e S\So are a consequence of the choices in Step (2). 

Step 1. We will make ip p canonical for each p G So in turn. Initially, we set S' = Sq. 
Take qeS' such that A q < Ap for all p e S'\ {q} and set S' = S'\ {q}. 

Consider ip q = Olio 1 PTa ■ To make v? q canonical, we wish to have a,; < ej iq /i jq 
for all 0 < i < r q . We will do this iteratively for i = 0, 1,..., r q . 

Case 1. at < e^qfi q . In this case, we do nothing. 

Case 2. m ^ ej, q /j iq and w p (0i, q (0)) ^ w p (0i+i, q (0)) for all p G S'. We replace 

g. f. i ^ i 

each 07’ q ,,q with a single 0 i+i, q in ip q . 

For all p e So, the p-valuation of will not decrease. In fact, for p G S' this 
is deduced from the condition of this case. For p G So\S', this is a consequence of 
Lemmas 15.171 and 15.181 

Case 3. a z > e.i A fi A and Wp o (0 ijq (0)) > u) Po (0i+i, q (0)) for some p 0 6 S'. In 

e . f. 

this case, we cannot simply exchange 07’ q “ ,q for 0 i+i, q without lowering the po- 
valuation of ifiq. 

Instead, take [ G S' the prime ideal that is “closest” to q; that is, /(q, [) ^ i{ q,p) 
for all p e S', and in the case of equality Aj ^ X p . By Lemma 15.171 and the election 
of 1 we must have i( q, 1 ) > i = ?'(po, q). 

We remove 0“* q from ipq and insert 0“j into ipi. If * < £ := 

*([, q) then 0 ^ q = 0 jj and g has not changed (we have only redistributed its p- 
parts). If i = £, we must check that w p (0^i(0)) ^ w p (0^, q (0)) for all p G Sq. For 
p G S' this is a consequence of Lemma f5.191 

By the maximality of 1(1, q) amongst all prime ideals in S', the relative situation 
of t q , t[, t Po in the non-optimised tree is as indicated in Figured] 
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Also, by the remarks following (fl5l) . we have necessarily </>^ q = </>(p 0 , q); hence the 
node t' in Figure [9] corresponds to a type in the optimised tree and A^ q = A q ° = A q . 
In particular, for any p e S with i(p, q) = £ we cannot have J(p, q) > 1(1, q). 

Now consider p e S 0 \S'. If p / q and /(p, q) < 1(1, q), Lemma [5.191 is also appli¬ 
cable and yields iv p (<j>i,i(0)) > w p (<t>e,q(0))- If d(p, q) > 1 (1, q), we have necessarily 
z(p, q) > i, as we have just remarked above. This clearly implies «(p, 1) = l. 

In both cases, p = q or i(p, q) > i( l, q) = l, Proposition 12.31 shows that: 


Wp (4>e, q (0)) 
w p (<j>i,i(8)) 


Vi + \i. p 
ei • • • e^_i 
Vt + X\ 

d ■ ■ • e£_i 


Vt + Af jq 

ei • • -et~i ’ 
Vt + * l q 

ei ■ ■ ■ ee-i 


In the last equality, we used A q = Min {A q , Af}, by the choice of q. Since A^ q = A q , 
we get w p (4>i,q(0)) = Wp ((/>e,i(0)) as desired. 


We continue in this way until = 0. 


Step 2. Having completed Step 1, all p-parts of g are canonical, so if D(g) = 0, 
then g e Ok(S'), completing the process after considering adequate choices of the 
Okutsu approximations ip p = <pp for all p e S\Sq . 

Step 3. If there exists q G Sq with degy) q > n q , we choose 1 e So such that 
I( q,0> /(q, p) for all p e Sq\ {q} and Aj 1 > Ap in the case of equality. 

Having chosen q and 1, the transfer occurs as follows for each 0 < * < r q + 1: 

f <j>i t i if i < i = i(q, 1), 

0i ’ q \^l mi if O € = *(q,t). 

By Lemma 15.191 the p-valuation of the resulting <p will not decrease, except 
possibly in the case where p = q. However, since </? q = <j ) q , we may increase the 
q-valuation of ^> q by choosing a better Okutsu approximation </> q to compensate 
any decrease in value. 

After this step we remove q from So, since tp q = </> q and then return to Step 1. 

As Step 3 reduces the number of prime ideals in So, this process will clearly end 
after at most #S iterations. □ 


6 . Proof of Theorem 13.31 

In the interest of clarity, we will only prove Theorem 13.31 in the case I = Ol- 
The proof in the case of an arbitrary fractional ideal is almost identical. Further 
details on the required adaptations can be found in [T51 Ch. 6]. 
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6.1. Precomputation. Recall that our set S of prime ideals has a total ordering 
S = {qi, ..., q s } satisfying (0. Denote ns '■= 2pes n P; w s '■= ws,o L and consider 
the following intervals of S: 

[a, 6] := {q., : a ^ j < b} , 1 ^ a < & ^ s. 

Definition 6.1. An order preserving partition of S is a decomposition S = I\ u 
• • • u It of S into the disjoint union of intervals Ij = [ aj,bj] with increasing end 
points b\ < ■ ■ ■ < bt. 

Take extended families of numerators goj y ■ ■ ■ , Qn ,,/, of Okutsu Ij -bases, for all 
1 < j < t, where nj := nj y That is, each g i j j has degree i, belongs to Ok {Ij) and 
wi j (ft,/..) is maximal amongst all rnonic polynomials in 0„[a;] of degree i. 

Consider multi-indices i = (ii,... ,it) of degree degi = ii + ■ ■ ■ + it and monic 
polynomials ft := g iuh ■ ■ ■ g H .i t e 0[x\. 

We may consider the version of MaxMin presented in Algorithm [21 


Algorithm 2 MaxMinfS 1 = Ji u • ■ • u It] algorithm 

Input: An order preserving partition S = I\ u • • • u I t of S c: p, and extended 
families {ftft : 0 < i < ni j } of numerators of Okutsu Tj-bases for all 1 ^ j < t. 
Output: A family i 0 , ii,.. ., i„ s e N 4 of multi-indices of degree 0,1,..., ns, re¬ 
spectively. 

1: i 0 <- (0,...,0) gN‘ 

2: for k = 0 —* ns — 1 do 

3: j <- Min {1 s; * ^ t : w h (ft J = w s (ft J} 

4: ife+1 ifc + Uj 

5: end for 


Such a decomposition of MaxMin will be useful for the proof of Theorem 13.31 

Definition 6.2. For indices 1 < a < b < s, we say that I = [a, b] admits precom¬ 
putation if after natural identifications, the algorithm MaxMlnfS 1 ] has the same 
output as 

(19) MaxMin [S = {qi} u ■ • ■ u {q Q _i} u/u {qb+i} u ■ ■ ■ u {q s }], 

where we consider the output o/MaxMin[J] as an extended Okutsu I-basis. 

By “natural identifications” we mean that if the k-th output of MaxMlnfS 1 ] is 
ifc = (iqi, • • •, iq„), then the k-th output of the algoritm is: 

jfc = (*qi! • • • j *q a —Ij */j *qt+l) • • • j *q s ) j 

while the ij-th output o/MaxMin[7] is (i qa ,... ,i qt ). 

The next result is an immediate consequence of the definition. 

Lemma 6.3. Let S = I\ u • • • u It be an order preserving partition of S . If all 
intervals Ij admit precomputation, then MaxMinfS' = Ji u • • • u It] has the same 
output as MaxMin [S], after natural identifications. □ 

Let us give a criterion for an interval to admit precomputation. 
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Lemma 6.4. Let io, ii,..., i ns be the output of MaxMinfS 1 ] and let I c: S be an 
interval of S. For each 0 < k ^ ns, let i& = (*q) qe g and denote 

9i k = | 1 9iq,qi Gi k = | | 9i q ,<\- 

qeS qeS\I 

Suppose that for each 0 ^ k < ns the following condition holds 

Wi (, g ik ) = ws (g ik ) => Wp (G ifc ) = w q (GiJ , V p,q e I. 

Then, I admits precomputation. 

Proof. Let (ife)o^fc^n s be the output of MaxMinfS] and (jfc)o^fe^n s be the output 
of the precomputed MaxMin algorithm (usd. 

Clearly, io and jo may be identified. For k ^ 0, suppose that i*, may be identified 
with jfc. This means 

ifc = (®qi ) • • • > *q s ) ) 

jfc = (*qi, • ■ • , iq a —1, */, *q b +l, • • • , iq 3 ) , 

while the ii -th output of MaxMin[/] is the multi-index (*q 3 -) a< . <6 - 

Let g ' 0 , ..., g' ni be the numerators deduced from the application of MaxMin[/] 
and g'- g , ..., gj the numerators deduced from (fl9l) . Clearly, 

b 

9jk = 9ij ]^[ gi q ,q = ]^5i,,q ]^[ <7*,,q = 9i k , 9i, = ]^[ 5i qm ,q m - 

q eS\I q el qeS\I rn=a 

The algorithm MaxMinfS] outputs i^+i = ifc + Uj, where 

3 = Min{l < m < s : w qm ( g ik ) = w s (g ik )} ■ 

If q j f /, then the q^-index in will also be the least index satisfying w qj (g'j k ) = 
ws (gj fc ), since g'- k = g, k . Thus, the algorithm in (fl9l> will also increase the q^- 
coordinate. 

If q j 6 /, then wi (g ik ) = ws (ffij and w qm (g ik ) > ws (gi k ) for all m < a; 
thus, (IT5T) will increase ij by one. In this case, we must show that the (zj + l)-th 
output of MaxMin[/] is the multi-index obtained from (iq,) a<7< /, by increasing the 
q.,-coordinate by one. 

The index increased by MaxMinf/] will be: 

J = Min {a < m < b : w qm (g-J = i uj (g- 7 )} . 

By hypothesis, v := w q (Gi k ) is independent of the choice of q e /. Since 
9i k = G ik g' n , we have: 

(ffij = wq (g' ; ) + v, Vq el. 

In particular, ws (gtj = Wi (g ik ) = wi (g'J + v, so that J = j. □ 

One specific case of precomputation which we will make use of, is the precom¬ 
putation of certain intervals S t S defined by a type t. 

Lemma 6.5. For any tsT, if the interval St is non-empty, it admits precompu¬ 
tation. 
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Proof. For every p e St and every q £ S t , the explicit formulas from Proposition 
12.31 show that w p is independent of p, for all i. Hence, the same is true for all 

polynomials G i k that are a product of these ^-polynomials. 

Thus, meets the criterion of Lemma 16.41 □ 

In [13] we give concrete examples of intervals / cz S which do not admitprecom- 
putation. 

6.2. The block-wise MaxMin algorithm. Consider an ordered subset S = 
{qi, ...,q,}£P. Let T s be the tree gathering all the paths of all leaves 4 P e T for 
p e S. Take 

£ = i(S) := Min {*(p, q) : p, q e S}. 

The tree Ts is disconnected if and only if £ = 0. Assume from now on that £ ^ 1; 
in this case, ^ — 1 is the order of the greatest common node of all paths joining 
the leaves of Tg with the root node. In this case, the Okutsu frames of all primes 
p G S have the same first £ — 1 key polynomials <f>i, ..., ■ Thus, the first me 

numerators of the Okutsu p-bases coincide for all p 6 S. Let 

J\f = {1 = ho, h ±,..., h m ^— i}, 
be the family of these common numerators. Note that 
(20) w p ( h ) = w q (h) , V p, q e S, V h e Af. 

Lemma 6.6. For all p, q e S and all 0 < r, t < m 

,Wq(h r+t ), if r + t < me, 


Wq (h r h t ) < 


e,phk), if r + t = mi + k, k ^ 0. 


Proof. By Lemma 15.171 (1). in any product of powers of <f>i,..., <f>e-i we may replace 


jei-ifi- 


□ 


H'i-i with (j>^p to increase the q-valuation. This proves both inequalities. 

Lemma 6.7. Let i be a maximal multi-index of degree divisible by mp. 

(1) There exists a maximal multi-index i' = (*p) p <=s of the same degree, having 
all its coordinates divisible by mg. 

(2) All elements in the family are maximal numerators. 

Proof. For 0 < j < me, let j = (j p ) p6 s be a multi-index of degree ime + j. Each 
index j p may be written 

j p = q p me + k p , 0 ^ fc p < me, 

and the numerators g ]p , v of the Okutsu p-basis may be written 

9j p ,p = G p h kf , degGp = q p m e . 

Since all polynomials G p have a degree which is a multiple of me, we have 
Spes = + j> f° r some non-negative integer q. 

Let g = n peS G p and choose any fixed prime ideal po G S. By (12U1) . an iterative 
application of the inequalities in Lemma 16.61 shows that for any q G S' we have 


(.9j) = 



= Wq {g) + ■ 



< w q (g) + ' 


'ho h > 


)■ 


Since this holds for all q e S, we deduce that ws (gj) < ws 


Iqhj 


y 
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These arguments, applied to j = i (and j = 0) prove item (1). Also, applied to 
an arbitrary j of degree deg(i) + j show that 

WS (ft) < ws (</>l Po ■ gj + w Po (. hj) < w s (ft) + Wp 0 (hj) = ws (fthj, 

by the maximality of ft. This proves item (2). □ 

Lemma 6.8. Let i = (i q ) qe s be an output multi-index of MaxMinfS 1 ] of degree 
divisible by me,. 

(1) All coordinates i q are divisible by me- 

(2) Let j = Min{l < m < s : u> qm (ft) = ws (ft)}- Then, the next me iterations 
of MaxMinfS'] increase the coordinate . 

Proof. All coordinates of io are zero; hence divisible by me- Thus, it suffices to 
prove that any output multi-index i = (* q ) qe s whose coordinates are all divisible 
by me satisfies (2). 

Let j = Min {1 ^ to < s : Wq m (ft) = ws (ft)}- If i = ifc is the fc-th output multi¬ 
index of MaxMin[5], the algorithm selects ifc+i = ifc + Uj. Since i qj is a multiple 
of me, we have g ik+1 = ft fe /ii; hence, 

(ftfc+i) = Wq (ftj + w q (/n) 5* ws (ftj + Wq {hi) = W S (ft fc+1 ) , 

for all q e S. Thus, w q (ft fc+1 ) = ws (ft fc+1 ) if and only if w q (ftj = ws (ft J. Thus, 
the next iteration increases the q-,-coordinate again. By iterating this argument, 
we get ft fc+ = gi k h me - 1 - At this point, the q^-coordinate will be increased once 
more to yield i k+m f = ifc + meUj. □ 

This result shows that MaxMinfS 1 ] works by blocks of length me. Thus, we may 
consider Algorithm [3j where we agree that me = 1 if Ts is disconnected. Note that 
for me = 1, MaxMin[S;mJ coincides with MaxMinU']. 


Algorithm 3 MaxMinfS 1 ; me] algorithm 

Input: An ordered subset S = {qi,...,qj cr 'p and extended families 

(ft,q : 0 < i < n q } of numerators of Okutsu q-bases of each q e S. 

Output: A family io, i mt , 12 m e , ■ ■ ■, i„ s / m( of multi-indices with deg ifc = k, having 
all coordinates divisible by me- 

1: io ■*— (0,..., 0) 

2: for k = 0 —* {ns/me) — 1 do 

3: j <- Min jl i ^ s : w Pz (ft fcm J = w s (ft fcm J } 

f(fc-fcl)m^ 4 ifcm^ T meMj 

5: end for 


Theorem 6.9. The output multi-indices of MaxMhJS 1 ; me] are maximal amongst 
all multi-indices of the same degree with coordinates divisible by me- 

Theorem 13.31 follows from Theorem 16.91 In fact, by Lemma 16.71 all output multi¬ 
indices of MaxMinfS 1 ; me] will be maximal and by Lemma 16.81 these multi-indices 
coincide with the output multi-indices of degree divisible by me of MaxMin[S']. 

Finally, Lemma 16.81 shows how to derive all intermediary output multi-indices of 
MaxMiiJS'] and Lemma 16.71 shows that these multi-indices are maximal too. 
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6.3. Branching cases. The basic idea for the proof of Theorem 16.91 is to split S = 
U vjD (U for “up” and D for “down”) into the disjoint union of two intervals which 
admit precomputation and then analyse the behaviour of MaxMinfS' = U u D] for 
which the multi-indices have only two coordinates. 

Lemma 16.81 shows that the output multi-indices of MaxMinfS'; me] coincide with 
the output of an ordinary application of the 2-dimensional MaxMin applied to 
the precomputations MaxMin[C/; m^] and MaxMin[T>; me]. We shall denote this 
algorithm by MaxMin [S' = U D;m/]. 

We distinguish four cases according to the structure of the tree Ts: 

Case (A). The tree T $ is disconnected, composed oft connected trees with root 
nodes ifo,i> ■ • ■, fio.t- We take D to he the connected component of^s with root node 
ipo,t- 

For Tg connected, the proof of Theorem 16.91 makes use of the structure of the 
non-optimised tree with base type te~i which is the greatest common node in all 
paths joining the leaves of Xs with the root node. 

Let fie be the first representative of te-i which leads to branching. Thus, before 
constructing fie, the Montes algorithm may have constructed other representatives 
of te-i admitting unibranch refinements. 

Let A m i n be the least slope (in absolute size) occurring in the branching based on 
fie- Let ^min c 5 be the subset of all prime ideals derived from branches of slope 
Amin of fie- 

Case (B). There exists a branch t with slope A m in which suffered refinement. In 
this case, we take D = St to he the set of all prime ideals derived from this branch. 
Note that fie,p ^ fie V p e D, and that there may be other \ m \ n -branches. 

.O 

t, i Of. .. ' 

yjQ •-• • • • •-•.o • • • o:; 

. A m i n , 

o t 


Figure 10. Case|(B)| Tree X‘“ p with at least one refined A m m-branch. 


Case (C). None of the Xmi n -branches suffered refinement, and there are no other 
slopes. In other words, A^ iP = A m in and (f>e. v = fie for all p e S. 

In this case, we take D = S t , for any choice of a A m i ^-branch t. 

Case (D). None of the A m j n -branches suffered refinement, but there are other 
slopes. In other words, A^ iP = A m in and fie )V = fie, for all p e S^m £ S. 

In this case, we take D = S'min- 

In all cases, we may change the ordering of V so that D and U = S\D are 
intervals. 


6.3.1. Proof of Theorem, 1 6'. .91 in cases ~(A )\ (B) and \(C)\ Denote 


Vf + A n 


c : = 


ei • • • e^_i 

0, 


if Ts is connected, 
if ‘Is is disconnected. 
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V’O*' 



(^£— I; ^minj )) 


(itl?— 15 (4>£i A m in; Ip*)) 


(^T—15 ((p£i A m in; VO) — t 


Figure 11. Case|(C)| Tree T s with only unrefined A m i n -branches. 


Vo*- 



(t^_ i; (</>", A", ip")) 


(tT_i; (cp' e ,\',ip')) 
(^e—i'i ((p£-> A m i n , 'ip*)) 


((p£, A m i n , VO) 


> X m i a -branches 


Figure 12. Case |i I ) j] Tree Is with with unrefined A m i n -branches 
and other slopes. 

The explicit formulas from Proposition 12.31 show that 

(21) w p ((pm, q) = ( m/me)c = w q ((p m ,p) , V p e U, q e D, V m ^ t. 

On the other hand, all ideas and criteria about precomputation apply to the 
MaxMin algorithms restricted to all multi-indices whose coordinates are divisible 
by me- Hence, m shows that D and U = S\D meet the condition of Lemma 16.41 
and both intervals admit precomputation. 

Denote the respective output families of numerators of MaxMin \U; me] and 
MaxMin[D; me] by: 

1> ffl> • • ■ > 9nu/mt i 1) 9\i ■ ■ • ) 9nn/mt :' 

Note that degg^ = kmi, for all k < njj/mi and deg g' k = kme , for all k < nn/me- 
By Lemma 15751 MaxMin[5; mi] has the same output as MaxMin[,S' = UkjD; me], 
after natural identifications of the respective multi-indices. In other words, if (i,j) 
is the fc-th output of MaxMinfS' = U D;me] (so that k = i + j), then the k- th 

numerator provided by MaxMin [S'; me] is g^g'. 

Definition 6.10. We say that a monic polynomial G e G[x] has support in a 
subset S' c S if it is a product of polynomials <p m ,p for p e S' and m ^ L 
Note that the degree of G is necessarily a multiple of me. 

In order to prove Theorem 16.91 we must show that the output numerators of 
MaxMinfS 1 ; me] are maximal amongst all polynomials of the same degree with sup¬ 
port in S. 

We proceed by induction on ffS. The case ffS = 1 being trivial, we may 
assume by the induction hypothesis that both sequences of numerators are maximal 
amongst all polynomials of the same degree with support in U and D , respectively. 
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For all 0 < z < njj/mt and all 0 ^ j < riD/me, denote 
( 22 ) Vi := wu (ffi) - ic, z/ := ( 3 ') - jc. 

We agree that z/_i = i /_ 1 = —1. 

Lemma 6.11. For all i,j ^ 0, 

t'i < *A+i, t'j < ^+i- 

Proof. By Proposition 12.31 w q (<j>e,q(8)) = (Vi + A^ q )/(ei • • • e^_i) for all q e U. 
Since A^ q ^ A m i n , the maximality of g l+ i implies 

wu (9i+ 1 ) ^ w v (gi<j>e,q) > % (&) + c, V q e 17. 

Similarly, the maximality of 3' +1 implies 

(Pj-+i) > w D (g'j&e) = w D ( 3 ') + c. 

By the definition (I221l of o,;, z/, this ends the proof of the lemma. □ 

For any bi-index i = (z, j), and any p G [/, q £ D, (12T1) shows that 
wu (3i) = wu ( 3i ) + jc = z/j + (z + j)c = Vi + (degi)c, 
wd (3i) = wd (Sj) + *c = i/'- + (z + j)c = v'j + (degi)c, 

WS (ffi) = W S (gig'j) = Min{z/j,z/} + (degi)c. 

Therefore, these numbers z/j, 1 / determine the flow of MaxMinfS 1 = U u D; me\. 
If (z, j) is an output pair, the next output pair is decided as follows: 

Wu ( gig'j ) = ws ( gig'j ) za < ij, “[/-minimal”, 

wd (, 9i9'j ) = w s ( 3 zSj) ^ “D- minimal”. 

The next output pair is (z + 1, j) in the [/-minimal case, and (i, j + 1) in the 
D-minimal case. 

Proposition 6.12. The output bi-indices ( i,j ) of MaxMin[5' = U u D\mi] satisfy 
the following properties: 

(1) Either vl_ 1 ^ z/, ; < l/j, or z/;_ 1 < v'j < 1 \. 

(2) The output multi-indices o/MaxMinfS 1 = U kj D;me\ are maximal amongst 
all polynomials of the same degree with suppose in S. 

Proof. Clearly, the initial output pair (0,0) satisfies (1). Let us check that if an 
output pair (i,j) satisfies ( 1 ), then the next output pair satisfies ( 1 ) as well. 


Suppose that z/'._ 1 

=% Vi 

^v', 

so that 

the next output pair is 

(* +1 , 3 ) 


v i+ i 


=> 

v 'o 

V/ 

1 

Vi 

/A 

+ 

/A 

Jn 



Vi+l 

>v 'o 

=> 

Vi 

/A 

JA 

< 

Vi+\. 


Suppose that za~i 

V 

V/ 

< Vi, 

so that 

the next output pair is 

(bj + !)• 


j: 

/A 

v 'j +1 

=> 

4 

< Vi 

< 

^+i> 



Vi > 

*4+1 

=> 

Vi- 

-1 < 

/ 

v i 

v' j+1 < Vi. 



This proves item (1). As a consequence, for any k 6 Z such that 0 ^ i — k < 
nu/mn and 0 < j + k ^ no/mi, we have: 

(23) Min { Vi_ k , v' j+k } < Min {v h i/j } . 
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In fact, if u'_ 1 < Vi < i/, then Min | u' +fe | < whereas in the case 

i'i-i Vj < Vi, we have Min j !/»_*,, u' +fc j < u'. 

In order to prove (2), suppose that (i, j) is an output pair of MaxMinfS” = 
U u D\me] and let g be a polynomial of degree (i + j)mt with support in S. We 
may write g = GG' , with G, G' polynomials with support in U and D , respectively. 

Suppose degG = (i — k)mt , degG' = (j + k)mt, for certain k e Z. By (12T1) and 
the maximality of the numerators gi-k, g'j+ki we have: 

wu ( g ) = wu (GG') = wu (G) + (j + k)c 

Ui_ fc + (i - k)c + (j + k)c = Vi- k + (i + j)c, 


w D {g) = w D (GG') = w D (G') + (i - k)c 

< Z'j+fc + (j + fc)c + (i - k)c = v' 3+k + (i + j)c. 
Hence, by using ([23]) . we get: 

ws (g) = Min {w v (g) , w D (g)} = Min v' j+k } + (i + j)c 
< Min {vi, v'j} + (i + j)c = w s (gig)). 

This ends the proof of Theorem 16.91 in cases |( A) (B) and |(C)| 


□ 


6.3.2. Precomputation in Case \(D j\ Recall that D = S m i n and U = S\D. In this 
case, we have: 

= <f>i = 0(p,q), V p e U,q e D. 

For each pe5we denote by A p the slope of the branch of in the non-optimised 
tree to which the leaf of p belongs. Also, we denote 


c := 


Vi + A n 


^ _ A p A n 

o n :— - 


Ss 0. 


The explicit formulas presented in Proposition 12.31 show that for all p e U, q e D: 
w q (0i, p ) = ( mi/me)c , 


(24) j (mi/m t )c, if i > £, 

W P W,q) = , .f ■ „ 

I Pp + C, it 1 = 1 . 

Let G be a polynomial of degree irng with support in U , and let G' be a polyno¬ 
mial of degree jrrif with support in D. If m := ord<^(G'), the formulas (1M1) show 
that: 

wu (GG') = Min {w p (G) + mS p } peU + jc, 
w D (GG') = w D (G') + ic. 

The first formula of pi shows that D meets the criterion of Lemma m and 
admits precomputation. In order to show that U admits precomputation too, we 
need another lemma. 


(25) 


Notation. For each pel?, we denote m p := p = Ci )P f( ]P mf. 

Lemma 6.13. Let i = (i p ) pe s be an output of MaxMinfS 1 ; me\ and g = gi the 
corresponding numerator. Let p e S be the least prime with w p (g) = ws ( g ). 
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(1) If p e D and mp \ ip, then the next ee,pfe,p output numerators are 


i> e / ,pfe,p \ 


and finally 


j^[siq,q I ' gip+m p ,p- 
\q^p / 


(2) If p e U, then m q | i q for all q e D. 

Proof. Suppose p e D and m p | ip. Since the element g* p , a numerator of the 
Okutsu p-basis, has degree divisible by m p , it is not divisible by <j>e, and c/j„+i, P = 
gi p ,p4>e- Hence, the next output numerator is g<fe- 

By (l24l) . Wp ( g<j>e ) = w p ( g) + c, while w q (g4>e) ^ w q (g) + c for all q e S. Thus, 
the least prime with w q (g<fe) = ws (g<f>e) is, once again, the prime p. 

This argument may be iterated as long as ord ^, e (gi p +km e ,p) = k < ee,pfe,p- For 
k = ee,pfe,p — 1 , the prime p is still the least one satisfying w p (g(f>g) = ws (gfie), 
so that the component of the multi-index corresponding to p is increased and the 
output multi-index is i + m p ui p . 

The second item follows immediately from the first. □ 

Corollary 6.14. U admits precomputation. 

Proof. Let us show that U meets the criterion of Lemma 16.41 

Let i = (ip)peS be an output of MaxMin[£; me] and let g = g^ be the corre¬ 
sponding numerator. Suppose that wjj (g) = ws (<?)■ With respect to the ordering 
of S, all elements in U are less than all elements in D; hence, the least prime p 
with Wp (g) = ws {g) belongs to U. By (2) of Lemma 16.131 m q | i q for all q e D, 
and this implies that none of the numerators 5 i q . q , for q e D, is divisible by eft,. 
Therefore, (l24ll shows that w v (<?i„, q ) = fi q /mt)c for all p e U, and the value 

Wp (Gi) = Wp (n q6C gi q , q ) is independent of p 6 U. □ 


6.3.3. Proof of Theorem \6.9\ in Case \(D~j\ Denote the respective output families of 
numerators of MaxMin[C/; me] and MaxMin[Zl; m<] by: 

1> 9li ■ ■ • ) gnu/mi 1 1> fll) • • • ) gno/me ' 

Note that deg < 7 *, = kmi , for all k < nu/me and deg g' k = km( for all k =% no/me. 

Let i = (i P ) P es be an output of MaxMlnfS 1 ; me]. Since U and D admit precom¬ 
putation, Lemma l6.3l shows that gi = gtg'j, for the fc-th output (i,j) of MaxMinfS 1 = 
U u D; me]. 

Notation. We denote [j] := ord fa(g'j), for 0 ^ j ^ no/me- 


By Lemma [6.131 all indices i q , for q £ D are divisible by m q except eventually 
for one, say i qo . Hence, [ 7 ] is the residue of the euclidian division of i qo by m qo . 
Note that [j] = 0 if and only if m q | i q for all q e D. 

Consider rational numbers z/,;, z/ as in 


|). The formulas (1251) translate into 
wu ( gig/ ) = Min {w p + [i ]^ P } pe(7 + jc, 

WD ( gig'j ) = W D (g'j) +ic = Uj + (i + j)c. 

Lemma 6.15. These data z/,;, v[. satisfy the following properties for all i,j > 0: 


(26) 
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(1) Vi-1 < Vi. 

(2) v'j_ 1 ^ v) and if [j] # 0 then equality holds. 

Proof. Take any pet/. By Proposition ^. 31 w p {4>e.,p) = (V^ + A^ p )/(ei ■ ■ • e^_i) > c, 
since A^ jq > A m i n . The maximality of gi implies 

wu {gt ) > w v {gi-i<j>i, p ) > wu {g^ i) + c. 

This proves (1). Similarly, the maximality of g) implies v)_j < v). 

By Lemma [6.131 [j] i=- 0 implies that g' :) = g)_i(j>t. Since w q (<pi) = c for all 
q efl, this implies Wd {g'j) = wd {g'j-i) + c - This proves (2). □ 

Lemma 6.16. Let i = {i,j) be an output pair o/MaxMin[5 = U u D: rile]. Then, 

(1) Either v)_ x ^ v t < v'j, or Vi-\ < v'j < z/j. 

(2) The next output pair is {i + 1, j) in the first case, and {i, j + 1) in the second 
case. 

Proof. Clearly, the initial pair (0,0) satisfies (1), the next output pair is (1,0), and 
it satisfies (1) too. Let us show by induction that if an output pair satisfies (1) 
then the next output pair is given as indicated in (2) and it satisfies (1) as well. 

Suppose that u'_ 1 < Vi ^ v). If the previous output pair was {i — 1, j), the 
induction hypothesis implies that we had v)_j ^ v,_ 1 =% v). Since v)_ ± ^ Vi_i < 
Vi < v'j , we have [j] = 0 by item (2) of Lemma |6.15l If the previous output pair 
was {i,j — 1), then z/'■_ 1 < Vi by the induction hypothesis. This leads again to 
v)_ 1 < v'j and to [)] = 0. Thus, (l26l) shows that 

wu {gig)) = wu {g^ + jc = Vi + {i + j)c 
< v ) + (* + j)c = W D {g l g)) . 

Thus, wu ( gig )) = ws {gig)) and the next output pair is (z + 1, j). The arguments 
of the proof of Proposition 16.121 show that {i + 1, j) satisfies (1). 

Suppose that z ^-i < v) < Vi. By (l26l) . we have 

WD {gig)) = v) + {i+ j)c <Vi + {i+ j)c 
= w v {gi) + jc = w v {gig)) ■ 

Thus, wd {did)) — w s {gig)) and the next output pair is {i,j+ 1). The arguments 
of the proof of Proposition 16.121 show that {i,j + 1) satisfies (1). □ 

Lemma 6.17. Consider indices 0 < k < i and let G be a polynomial of degree 
{i — k)mi with support in U. Then, 

Min {w p (G) + kd p } peU < Vi + {i - k)c. 

Proof. Let q e U be a prime ideal with a maximal value of A q . The statement 
follows from the following chain of inequalities: 

(27) Min {w p (G) + k5p} peU + kc wu {G(j^ q ) < wu {gi) = Vi + ic. 

The second inequality of (1271) follows from the maximality of gi. The first in¬ 
equality is deduced from the formulas from Proposition ^. 31 In fact, for any pel/, 
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these formulas yield wp {<f>e,q) = (Ve + A)/(ei • • • e^_i), for a certain rational number 
A, depending on p, such that A ^ A p ; hence, 

w p (G<j)f J = wp (G) + k—^—^— 
ei---ee~i 

^ Wp (G) + k(c + Sp), 

for all p 6 17, which implies the first inequality in (1271) . 

More precisely, if z(p,q) = £, then A = A p or A = Min{A p ,Aq}, according to 
</>(p, q) being equal to ^, q or not. Now, if p and q belong to the same ^-branch of 
the non-optimised tree, we have (see Definition 15.111) 

A p = A q < Min{A£, Ajj} < A. 

If p and q belong to different ^-branches of the non-optimised tree, then A p = A p , 
Aq = Aq, so that, again, 

A p = Min{A p , Aq} = MinjAj], A^} A. 

Finally, if i{ p, q) > £, then A = A^q = A^ p > A p , by Remark l5.12l □ 

We are ready to prove Theorem 16.91 in Case |(D)| 

Proposition 6.18. In Case [TtTJ} all output multi-indices of MaxMin},?; mi] are 
maximal amongst the multi-indices of the same degree whose coordinates are all 
divisible by me- 

Proof. Let i = (ip) pe s be an output multi-index of MaxMinfS 1 ; mg], obtained by 
joining the i-th output of MaxMin[[7] and the j-th output of MaxMin[77J. 

Let g be a polynomial of degree (i + j)rrif with support in S. We may write 
g = GG', with G, G' polynomials with support in U and D 1 respectively. 

Suppose degG = (i — k)mg, deg G' = (j + k)mg, for certain k e Z. Let us write 

G' = Hffi, fa f H, deg H = qm e . 

Note that q + m = j + k. By (USD, 

wu (GG') = Min {w p (G) + rnS v } peU + (j + k)c, 

w D {GG') = w D (G') + (* - k)c. 

Since w p {<f>i) = c for all q e D, the last equality leads to 

wd {GG') = wd (H) + (m + i — k)c < wd {g' q ) + (m + i — k)c 
= v' q + {q + m + i — k)c = v' q + (i + j)c. 

By Lemma 16.161 we may distinguish two cases according to the comparison of 
Vi with vl. 

Case 1. v'-_ 1 < Vi ^ z/. In this case, we saw during the proof of Lemma 16. 161 that 
[j] = 0. Hence, ws {di9j) = {didj) = v i + (* +i)c, by (1761) . We want to show 
that 


ws {GG') = Min {wu {GG') , w D {GG')} ^ v t + (i + j)c. 
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If m ^ k, then Lemma Th. 1 71 shows that 

wjj (GG') = Min {u> p (G) + m8 p } p€U + (j + k)c 

< Min {w p (G) + k8 p } peU + (j + k)c 

< Vi + (z — k)c + (j + k)c = Vi + (z + j)c. 

If to > k, then q < j , or equivalently q ^ j — 1. Thus, (l28l) shows that 
w D (GG') < v' q + (i + j)c < + (i + j)c < Vi + (i + j)c. 

Case 2 . Vi_, <v'< Vi. In this case, we saw during the proof of Lemma 16.161 that 
W S (di9j) = w d (, g%gj ) = v'j + (i + j)c. We want to show that 

w s (GG') = Min{ Wy (GG') ,w D (GG')} < v) + (i + j)c. 

If to < k , then to ^ fc—1. Having in mind that deg G/mi = i—k = (i — 1) — (k — 1), 
Lemma T6 .1 71 shows that 

wu (GG') = Min {w p (G(6)) + m5 p } peU + (j + k)c 

< Min {w p ( G(6 )) + (k- l)^ P } peC7 + (j + k)c 
^ i + (i — fc)c + (j + fc)c = v ^i + (z + j)c ^ u' + (z + j)c. 

If to ^ fc, then g ^ j and (1251) shows that 

w D (GG') < u' + (z + j)c < u' + (z + j)c. 

□ 
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